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The recursion method is used to determine the T = spin dynamic structure factor S zz (q,u>) 
in the Luttinger liquid state and in the phase-separated state of the one- dimensional t-J model. 
As the exchange coupling increases from zero, the dispersions and line shapes of the dominant spin 
excitations are observed to undergo a major metamorphosis between the free-fermion limit and the 
onset of phase separation. The familiar two-spinon spectrum of the Heisenberg antiferromagnetic 
chain emerges gradually in the strongly phase-separated state. 
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The t-J model for strongly correlated electrons has 
been widely used as a prototypical system for the study 
of the iaterplay between charge and spin degrees of 
freedomJd For a one-dimensional (ID) lattice of N sites 
the Hamiltonian reads 



H t -j = -t + 4fi,<A<r} 

N r ii 

+ JJ2 \Si-S l+1 --n ini+ A (1) 
1=1 ^ - 1 

with ci,a = Q )CT (l-n;,_ CT ), ni = m^ + ni^, n/ )ff = c\ a c^ a , 

&l = \{ n iA ~ n i-l)> an d S/~ = c[^ci,i- Here we consider 
the quarter- filled-band case (N e — N/2 electrons). 

The hopping term amounts to an effectively repulsive 
force between electrons on nearest-neighbor sites,0 and 
the exchange term represents an attractive force if these 
two electrons have opposite spins. As the exchange in- 
teraction increases from zero, the ground state of H t -j 
undergoes a transition, at J ft ~ 3.2, from a Luttinger liq- 
uid state to a phase-separated state.B~El In the Luttinger 
liquid state, both the charge correlations (at q = ir) and 
the spin correlations (at q — ir/2) are critical. The tran- 
sition, which is driven by the spin coupling, produces 
charge long-range order (at q = 0) combined with a new 
type of spin criticality (at q — n). 

The focus of this study is on the frequency-dependent 
spin fluctuations of H t -j as they manifest themselves in 
the T — spin dynamic structure factor 
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For the calculation of this quantity, we employ the re- 
cursion methods in combination with a strong- coupling 
continued-fraction analysis.uB The recursion method in 
the present context is based on an orthogonal expansion 
of the wave function | = S*(-t)\G), where S z q = 

spin fluctuation operator, and \G) 
is the finite-size ground-state wave function of (|l|). Af- 
ter some intermediate steps, the algorithm produces a se- 
quence of continued- fraction coefficients Af (q), A^q), . . . 
for the relaxation function, 
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which is the Laplace transform of the symmetrized cor- 
relation function ft(S*(t)S*)/(S*S^ q ). The T = spin 
dynamic structure factor (g) is then obtained from (||) 
via the relation 

S zz (q, u) = i(S*S z _ g )Q(uj) lim »[cg* {q,s - iu)]. (4) 



In a previous paperi we have used the recursion 
method together with a weak- coupling continued-fraction 
analysis for a study of the charge dynamic structure fac- 
tor in the Luttinger liquid phase of the t-J model. The 
spin dynamics poses a far greater challenge. Very few 
explicit results seem to exist.EEl One key spin dynamical 
property in the Luttinger liquid phase can be inferred 
from the asymptotic behavior of the-static spin correla- 



tions as proposed in previous work! 
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where rj p is the exponent which also governs the alge- 
braic decay of the (q — 7r)-oscillations, ~ cos(7^m)/m' ,f, , 
in the static charge correlation function (nini +rn ). This 
exponent is known to assume the value r\ p — 2 in the 
free-fermion limit (J/t = + ). It increases, linearly, 
An p ~ 0.40 J/t, in the weak-coupling regime,Ha assumes 
the value rj p ~ 3.4 at J/t = 2 (supersymmetric case), Ell 
and diverges,at the endpoint, J/t ~ 3.2, of the Luttinger 
liquid state.0 

The oscillatory term in (|^) implies that the dynami- 
cally relevant excitation spectrum of S zz (q,w) is gapless 
at q = 7r/2. The spectral-weight distribution of the spin 
dynamic structure factor at this critical wave number 
can then be predicted (under mild assumptions) to have 
a singularity of the form 



^(tt/2,^)^^/ 4 - 1 
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The infrared exponent starts out negative in the free- 
fermion limit, ~ w -1 / 2 , increases monotonically with in- 
creasing J/t, and then diverges at J/t ~ 3.2, where phase 
separation sets in. A landmark change in S' zz (7r, u)) is ex- 
pected to occur at the point where the infrared exponent 
switches sign (from negative to positive). This happens 
for rj p = 4, which corresponds to the coupling strength 
J/t ~ 2.3. 

Our results for S zz (q,uo) indicate that the Luttinger 
liquid phase of the t-J model can be divided into two 
regimes with distinct spin dynamical properties. As a 
representative result of the first regime (0 < J/t < 1), 
we show in Fig. 1 S zz (q,uj) as a continuous function of 
lj and a discrete function of the g-values realized in a 
system of N — 12 sites with coupling strength J/t = 0.1. 

All results presented here have been calculated via 
a strong-coupling continued-fraction reconstruction from 
the coefficients Ai, . . . , Ag and a Gaussian terminator. 
The Afc's have been extracted via the recursion method 
from the ground-state wave function for N = 12, which 
in turn has been computed via the conjugate-gradient 
method. The entire procedure was explained in Ref. |?]. 

Throughout the Brillouin zone except at small q we 
observe a well-defined dynamically relevant spin mode 
with a | cos q \ -like dispersion as indicated by the full 
circles. Near the critical wave number, q = n/2, the 
function S zz (q,Lo) may exhibit a power-law divergence 
of the form - [uj 2 - c 2 {n/2 - g) 2 ]^/ 4 - 1 )/ 2 , similar to 
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what has been observeeLat the critical wave number in 
other Luttinger liquids £3 At long wavelengths the spec- 
tral weight in S zz (q,uj) is concentrated at fairly low fre- 
quencies. Data for longer chains will be needed for a 
quantitative analysis of the spin dispersions at small q in 
this regime. 

As the exchange coupling increases toward J/t ~ 1, 
the following changes in the spectral-weight distribution 
of S zz (q,Ld) can be identified: 

(i) The amplitude of the | cos g|-like dispersion grows 
with increasing J/t. The gradual upward shift of the 
peak position in S zz (q, u) is accompanied by a significant 
increase in linewidth. For the q — -k spin mode this 
is contrary to what one expects under the influence of 
an antiferromagnetic exchange interaction of increasing 
strength. That trend changes at stronger coupling as we 
shall see. 

(ii) The intensity of the central peak in S zz (tt/2,lo) 
weakens gradually. This is illustrated in the main dia- 
gram of Fig. 2. The peak turns shallow and then vanishes 
altogether. This property of our data reflects the grad- 
ual weakening and ultimate disappearance of the power- 
law divergence (Jfy , given the (approximately known) J/t- 
dependence of the infrared exponent. 

As the coupling strength increases past the value J/t ~ 
1, the spin modes which dominate S zz (q,u>) in the first 
regime of the Luttinger liquid phase broaden rapidly and 
lose their distinctiveness. There is a crossover region be- 
tween the first and second regimes, which roughly com- 
prises the coupling range 1 < J/t < 2. Over that range, 
the spectral weight in S zz (q, to) is distributed over a broad 
structure with rapidly shifting peaks. 

At the end of the crossover region, a new type of spin 
mode with an entirely different kind of dispersion has 
gained prominence in S zz (q,u>), and it stays dominant 
throughout the remainder of the Luttinger liquid phase. 
This is illustrated in Fig. 3 by a plot of S zz (q,w) for a 
J/i-value near the onset of phase separation. The repre- 
sentation is the same as in Fig. 1 except for the different 
frequency scale. 

The dispersion of these new spin modes, as shown 
by the full circles, does no longer have a soft mode at 
q = tt/2. At J/t ~ 2.0 the dispersion has a |sin(q/2)|- 
like shape, i.e. a smooth maximum at q = -k and a ten- 
dency to approach zero linearly as q — > 0. As J/t in- 
creases toward the transition point, the peak positions in 
S zz (q,uj) gradually shift to lower values of u/J (with t 
held fixed). The shift proceeds more rapidly at q near 
7r than at smaller q, which has the consequence that 
the maximum in the dynamically relevant spin disper- 
sion starts to move away from q = tt toward q = n/2. 
This trend from a | sin(q/2)|-like toward a | singj-like dis- 
persion persists in the phase-separated state (J/t > 3.2). 

The evolution of the spectral-weight distribution in 
S zz (jr,u}) from a well-defined mode at u/J ~ 1.4 near 
the onset of phase separation toward a soft mode in the 
strongly phase-separated state is shown in the inset to 
Fig. 2. By contrast, the spectral-weight in S zz (ir/2,u), 



which has been slowly shifting toward lower values of 
lu/J in the second regime of the Luttinger liquid phase, 
now starts to move back out to higher frequencies in the 
phase-separated state (not shown here). The single-peak 
structure grows taller, the linewidth shrinks somewhat, 
and the peak position settles at to/ J ~ 1.57. 

All these properties and trends reflect or foreshadow 
the much better understood T = dynamical prop- 
erties of the ID s = 1/2 Heisenberg antiferromagnet 
- the end product of the t-J model in the completely 
phase-separated state. The spin dynamic structure fac- 
tor S zz {q,w) under these circumstances is known to be 
dominated by a continuum of two-spinon excitations 
with a lower boundary £l(<?) = (tt/2)J\ sinq|, where the 
spectral-weight distribution has a divergent singularity, 
and an upper boundary eu(q) = 7rJ| sin(g/2)|, whem^the 
spectral-weight distribution becomes very smalliJfJEj 
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FIG. 1. Spin dynamic structure factor S zz (q,co) at T = 
for q = 27T//12, 1 = 0, . . . , 6 of the ID t-J model with t = 1 and 
J = 0.1 in the Luttinger liquid phase near the free-fermion 
limit. 

FIG. 2. Line shape at q = tv/2 of the spin dynamic struc- 
ture factor S zz (q,u>) at T = for various values of J in the 
Luttinger liquid state of the ID t-J model with t = 1 (main 
plot). Line shape at q = ir for various values of J in the 
phase-separated state (inset). 

FIG. 3. Spin dynamic structure factor S zz (q,io) at T — 
for g = 27TZ/12, Z = 0, . . . , 6 of the ID t-J model with t = 1 
and J = 2.75 in the Luttinger liquid phase near the onset of 
phase separation. 
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S zz (q,cj) 




JS 77 (q,cj) 




S 77 (tt/2,cj) 




